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ABSTRACT 

A  generalized  one  dimensional  computer  code 
for  analyzing  the  flow  and  heat  transfer  in  the  tur¬ 
bomachinery  cooling  passages  has  been  developed. 
This  code  is  capable  of  handling  rotating  cooling 
passages  with  turbulators  (inline,  staggered,  or  in¬ 
clined),  180  degree  turns  (both  sharp  and  gradual 
turns,  with  and  without  turoulators),  pin  fins,  finned 
passages,  by-pass  flows,  tip  cap  impingement  flows, 
and  flow  branching.  The  code  is  an  extension  of  a 
one-dimensional  code  developed  by  P.  Meitner.  In 
the  subject  code,  correlations  for  both  heat  transfer 
coefficient  and  pressure  loss  computations  have  been 
developed  to  model  each  of  the  above-mentioned  type 
of  coolant  passages.  The  code  has  the  capability  of 
independently  computing  the  friction  factor  and  heat 
transfer  coefficient  on  each  side  of  a  rectangular  pas¬ 
sage.  Either  the  mass  flow  at  inlet  to  the  channel  or 
the  exit  plane  pressure  can  be  specified.  For  a  spec¬ 
ified  inlet  total  temperature,  inlet  total  pressure  and 
exit  static  pressure,  the  code  computes  the  flow  rates 
through  the  main  branch  and  the  sub  branches  (if 
any),  flow  through  tip  cap  for  impingement  cooling 
(if  any),  in  addition  to  computing  the  coolant  pres¬ 
sure,  temperature  and  heat  transfer  coefficient  distri¬ 
bution  in  each  coolant  flow  branch.  Predictions  from 
the  subject  code  for  both  non-rotating  and  rotating 
passages  agree  well  with  experimental  data.  The  code 
w.as  used  to  analyze  the  coolant  passage  of  a  research 
cooled  radial  rotor. 

INTRODUCTION 


higher  turbine  inlet  temperatures  are  desired.  As  the 
cooling  scheme  becomes  more  complex,  a  better  un¬ 
derstanding  of  the  internal  pressure  and  heat  transfer 
coefficient  distributions  is  necessary  to  be  able  to  pre¬ 
dict  the  hot  section  metal  temperatures  accurately. 
Typical  cooling  concepts  of  modern  turbine  blading 
are  shown  in  figure  1.  Figure  1(a)  shows  the  cooling 
scheme  for  an  axial  rotor,  while  figure  1(b)  shows  a 
hypothetical  cooled  radial  rotor.  The  shown  cooling 
passages  may  have  heat  transfer  enhancement  devices 
in  the  form  of  turbulators,  pin  fins,  finned  passages, 
leading  edge  tip  cap  impingement  flows,  etc.  The 
main  coolant  flow  may  also  be  split  into  sub  branches 
as  shown  in  figure  1(a)  and  1(c).  It  should  be  noted 
that  to  obtain  a  2-D  distribution  of  heat  transfer  co¬ 
efficient  in  the  plane  of  the  blade  surface  will  really 
require  a  3-D  analysis.  Because  of  the  complex  nature 
of  the  flow  field  in  cooling  passages,  it  is  extremely  dif¬ 
ficult  to  develop  a  generalized  3-D  flow  analysis  code 
capable  of  modeling  different  cooling  passage  config¬ 
urations  and  obtaining  detailed  distribution  of  heat 
transfer  coefficients  and  pressure  drops.  Such  infor¬ 
mation  is  required  in  the  detailed  thermal  and  struc¬ 
tural  analysis  of  the  rotor  and  blading  [2].  Hence  it 
was  felt  that  development  of  a  generalized  1-D  code 
which  can  model  various  types  of  cooling  passages  in¬ 
cluding  flow  branching  should  be  useful  in  the  design 
and  analysis  of  cooled  turbomachinery.  Based  on  the 
literature  review,  the  authors  conclude  that  such  a 
generalized  1-D  code  is  not  available  in  the  open  lit¬ 
erature. 


Increasing  the  turbine  inlet  temperature  is  a 
very  lucrative  approach  to  improving  the  thermal  ef¬ 
ficiency  and  specific  thrust  of  gas  turbine  engines. 
However,  this  practice  is  limited  by  the  availability 
of  materials  which  can  withstand  these  high  temper¬ 
atures.  Current  materials  cannot  be  exposed  to  tem¬ 
peratures  above  1250  K  without  a  dramatic  decrease 
in  life  [1].  Blades  must  be  cooled  if  significantly 


The  code  reported  in  this  paper  is  an  extension 
of  the  code  developed  by  Meitner  [3],  and  is  appli¬ 
cable  to  any  combination  of  cooling  passage  config¬ 
urations  mentioned  above.  Specifically,  the  reported 
code  differs  from  the  code  of  reference  [3]  in  that  it 
can  calculate  and  balance  flows  in  sub-branches,  and 
can  thus  compute  flows  in  passages  shown  in  figure 


1(a)  and  1(c)  wliereag  tlie  code  of  reference  [3]  will 
not  be  capaole  of  handling  these  passages.  Also, 
the  code  reported  in  this  paper  can  make  separate 
calculations  for  the  lcading(suction)  surface,  trail- 
ing(pressure)  surface,  and  side  walls  of  a  rotating 
channel  of  rectangular-  cross  section.  This  is  very 
important  for  capturing  the  experimentally  observed 
behavior  in  rotating  environment.  In  addition,  a 
more  detailed  modeling  of  each  type  of  cooling  pas¬ 
sage  above  has  been  incorporated  in  the  present  code 
compared  to  the  code  of  Meitner.  The  code  can  pre¬ 
dict  the  mass  flow  rates,  internal  pressure  and  tem¬ 
perature  distribution,  heat  transfer  coefTicient  and 
lieat  flux  distribution  at  the  walls  along  the  passage. 
In  a  limited  sense,  it  can  predict  the  heat  transfer 
coeflicient  distribution  along  each  side  of  rectangu¬ 
lar  passage.  It  integrates  the  1-D  compressible  mo¬ 
mentum  and  energy  equations  along  a  defined  flow 
path  and  will  account  for  area  changes,  mass  ad¬ 
dition/subtraction,  Coriolis  acceleration,  centrifugal 
force  effects,  bouyancy  force  effects,  friction  and  heat 
addition.  Flow  can  be  bled  off  for  tip  cap  impinge¬ 
ment  cooling  and  a  flow  ’bypass’  can  be  specified  in 
which  the  coolant  flow  is  taken  off  at  one  point  in 
the  flow  channel  and  re-introduced  at  a  point  fur¬ 
ther  downstream  in  the  same  channel.  As  mentioned 
above,  other  types  of  passages  that  can  be  modeled 
include  plain  passages,  rib-roughened  passages  (stag¬ 
gered  and  inline  turbulators  at  various  angles  of  at¬ 
tack),  dilferent  pin  fin  arrangements  (with  or  without 
lateral  flow  ejection,  as  shown  in  figure  1(a)),  finned 
passages,  180  degree  turns  (both  sharp  and  gradual) 
and  flow  brandling. 

n.ased  on  the  most  recent  experimental  data 
available,  correlations  have  been  obtained  to  model 
the  various  types  of  rotating  cooling  passages.  The 
user  has  the  option  of  specifying  the  temperature  and 
pressure  at  inlet  to  the  main  branch  and  either  the 
main  flow  rale  or  the  exit  pressures  from  the  sub 
branches.  The  present  version  of  the  code  cannot 
accomodate  su[>ersonic  flow  anywhere  in  the  cooling 
pa.ssage.  The  flow  can  at  best  be  choked  at  exit  from 
one  or  more  of  the  branches.  However,  it  is  felt  that 
this  is  not  a  serious  restriction  for  modeling  turboiiia- 
chincry  cooling  passages. 

Comparisons  have  been  made  with  a  wide  vari¬ 
ety  of  non-rotating  and  rotating  passage  experimen¬ 
tal  data  representing  various  types  and  combinations 
of  cooling  passages.  It  is  shown  that  the  predictions 
compare  well  with  the  experimental  data. 


(a):  Axial  Turbine  [9] 


(b):  Radial  Turbine 


(c):  Radial  tubine  with  Flow  branching 

Figure  1.  Typical  ('ooling  Concepts  of  Modern  Gas 
Turbines 
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FLOW  PATH  MODELING  AND 
COMPUTATION  PROCEDURE 

The  flow  path  is  defined  in  cartesian  co-ordinates 
i,r  and  y  as  shown  in  figure  1(b).  The  flow  path  ge¬ 
ometry  is  described  at  nodes  and  the  type  of  passage 
between  each  set  of  nodes  is  specified  using  different 
parameters  in  the  code  input. 

For  specified  entrance  conditions  and  a  specified 
(or  guessed)  flow  rate,  the  flow  solution  is  obtained 
by  marching  along  the  defined  flow  path  from  node 
to  node.  The  pressure  and  temperature  at  each  node 
is  calculated  from  the  values  at  the  previous  node  by; 


1  /  thRT  \  \  dA 


CpA\  pA 


1  /  mRT 


pCp I  pA 


4hcA(T^-T^^) 

fnCJp  Dfi 


where 


Pn  =  Pn-1  +  {dp/dx)n  *  dXn 


(dm/dx)n  =  (mn  -  m„_i)/(Z„) 


r„  =  r„_i  +  (dr/dx)„*dx„  (2) 

where  (dp/dx)n  and  (dT/dx)n  are  given  by  the  solu¬ 
tion  of  momentum  and  energy  equations  as  [3]: 
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The  constants  Cl,  C2  and  CT  are  defined  as  fol¬ 
lows  [3]: 
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Iteration  takes  place  in  each  marching  step  until 
a  convergence  is  acheived.  After  marching  through  all 
the  nodes,  the  calculated  branch  or  sub  branch  exit 

firessures  are  compared  with  the  given  exit  pressures 
if  specified).  If  necessary,  a  correction  is  made  to  the 
inlet  flow  rate  and  the  marching  procedure  repeated 
until  the  calculated  and  specified  exit  pressure  values 
agree  within  a  prescribed  tolerance. 

For  flows  where  the  local  coolant  passage  Mach 
number  becomes  high  (>  0.9),  pressure  changes  in  the 
flow  direction  become  highly  non-linear.  A  provision 
has  been  made  in  the  code  to  automatically  add  ad¬ 
ditional  nodes  between  the  original  nodes  and  repeat 
the  iteration  process  mentioned  above  until  the  Mach 
number  values  at  each  of  the  original  nodes  do  not 
change  between  successive  increases  in  the  number  of 
nodes. 

In  the  present  paper,  for  the  different  types  of 
cooling  passsages  mentioned  above,  both  friction  fac¬ 
tor  and  heat  transfer  coefficients  are  modeled  with 
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correlations  developed  from  experimental  data  avail- 
able  in  the  literature  or  using  the  correlations  avail¬ 
able  in  the  literature  (with  modifications).  For  val¬ 
idating  the  correlations  developed,  computations  us¬ 
ing  these  correlations  have  been  compared  with  ex¬ 
perimental  data  from  an  entirely  different  source  than 
those  used  in  the  development  of  correlations. 

The  correlations  developed  are  generally  of  the 
form: 


Nu  or  f  ■a(NDiy‘{ND2Y{ND3y...  (11) 

where  Nu  is  the  Nusselt  number,  f  is  the  friction 
factor,  NDl,  ND2,  etc  are  relevant  non-dimensional 
quantities  like  Reynolds  number  (Re),  Rotational 
Rayleigh  number  (Ra),  Rossby  number  (Ro),  Prandtl 
number  (Pr),  rib  licight  (or  roughness  heignt)  to  hy¬ 
draulic  diameter  ratio  (a/Dh),  the  rib  pitch-to-height 
ratio  (P/e),  channel  aspect  ratio  (W/H),  pin-fin  span- 
wise  spacing  to  diameter  ratio  (s/d),  pin-fin  stream- 
wise  spacing  to  diameter  ratio  (X/d),  pin  height  to 
pin  diameter  ratio  (H/d),  etc.  The  correlations  were 
obtained  from  the  experimental  data  using  the  SAS 
[4]  multiple  correlation  and  regression  analysis  pro¬ 
cedure.  The  correlation  equations  are  not  shown  in 
the  paper  due  to  space  limitations.  Additional  de¬ 
tails  of  the  correlation  equations  and  the  code  will  be 
presented  in  a  future  NASA  report. 

In  the  present  version  of  the  1-D  code,  correla¬ 
tions  have  been  included  to  model  any  combination 
of  the  following  types  of  passages: 

a)  Plain  pa.ss.age3  (Smooth  and  Rough):  Based  on 
the  magnitude  of  relative  roughness  parameter 
(c/D/,),  the  passages  are  classified  as  smooth  or  rough 
and  the  corresponding  correlation  equations  are  used. 

I'Viction  Factor;  For  laminar  flow  in  circular,  rectan¬ 
gular  or  square  ducts,  based  on  the  data  available  in 
Kays  and  Crawford  [5]  in  the  form  of  tables  and  fig¬ 
ures,  correlations  have  been  obtained  to  model  both 
smooth  aTul  rough  non-rotating  passages  and  for  both 
fully  developed  and  entrance  regions  of  flow.  For  tur¬ 
bulent  flow  and  in  the  transition  region  in  ducts,  cor¬ 
relations  listed  in  Welty,  Wicks  and  Wilson  [6)  have 
been  n.ned  for  both  smooth  and  rough  passages 

Heat  Transfer  CoefTicient:  For  laminar  flow,  correla¬ 
tions  are  developed  based  on  the  data  of  Kays  and 
Crawford  [5]  for  the  fully  developed  as  well  as  the 
entrance  regions  of  flows  in  both  circular  and  non- 
circular  ducts.  For  fully  developed  turbulent  flows, 
correlations  given  bv  Krcith  and  Bohn  [7]  and/or 
those  in  reference  (5]  liave  been  included.  Entrance 
effects  are  modeled  based  on  Kays  and  Crawford  [5]. 

b)  Turbulated  Passages:  Correlations  have  been  in¬ 
corporated  to  model  the  heat  transfer  coefficient  and 
pressure  drop  dl  tributions  for  rectangular  or  circular 
passages  with  transverse  rib  rougheners  (both  inline 
and  staggered  arrangements  and  at  any  given  angle 
of  attack).  Correlations  have  been  included  for  both 
the  turbulated  walls  and  the  smooth  side  walls  based 


on  the  work  of  Han,  et  al.  [8,9|  for  in  line  rib  ar¬ 
rangement  and  Taslim  and  Spring  [lO]  for  staggered 
a^angement.  Correlations  include  the  effects  of  non 
dimensional  parameters  like  rib  height  to  hydraulic  di¬ 
ameter  ratio  {efDh),  rib  pitch-to-height  ratio  (P/c), 
channel  aspect  ratio  (W/H).  Reynolds  number  (Re) 
and  rib  angle  of  attack  ratio  (a/90).  Based  on  Boyle’s 
data  [11],  entrance  effects  have  also  been  incorporated 
in  the  correlations. 

c)  Pin  Fins:  Based  on  the  work  of  Metzger,  et  al. 
[12,13,14,15],  correlations  arc  available  to  model  pres¬ 
sure  drop  and  heat  transfer  coefficients  (array  aver¬ 
aged)  for  different  staggered  pin  fin  configurations. 
Nusselt  number  (Nu)  correlations  which  include  the 
effects  of  streamwisc  pin  spacing  are  given  by  Metzger 

(12) .  The  local  row  resolved  heat  transfer  coefficient 
correction  has  been  correlated  based  on  the  data  of 
Metzger  [12]  as  well  as  the  corrections  to  Nu  to  in¬ 
clude  effects  of  channel  convergence.  Nusselt  number 
correlations  are  also  available  based  on  the  work  of 
VanFossen  [16].  It  should  be  noted  that  in  VanFos- 
sen  correlations,  the  Nusselt  numbers  are  based  on 
hydraulic  diametet  DV  which  depends  on  the  open 
volume  in  the  pin  fin  bank  and  the  total  heat  trans¬ 
fer  surface  area  whereas  the  Metzger  correlations  arc 
based  on  pin  diameter.  In  a  recent  review  paper  on 
staggered  array  pin  fin  heat  tansfer,  Armstrong  and 
Winstanley  [17]  recommend  the  use  of  Metzger  or 
VanFossen  correlations  for  short  pin  fins  (Il/d  <=  3) 
and  the  correlation  of  Faulkner  [18]  for  long  pin  fins 
(H/d  >  3).  Friction  factor  correlations  fi  r  both  short 
and  long  pin  fins  are  ba.scd  on  the  work  of  Metzger 

[13) .  ^ 

The  effect  of  channel  convergence  on  friction  fac¬ 
tor  is  incorporated  based  on  the  work  of  Brown  [19]. 
For  pin  fin  channels  with  lateral  ciection,  the  correla¬ 
tions  developed  by  Lau  et  al.  [20f  arc  used. 

d)  180  Degree  Turns:  Based  on  the  experimental  data 
of  Metzger,  et  al.  [21,22,23,24],  correlations  were  de¬ 
veloped  for  heat  transfer  coefficient  and  pressure  loss 
in  smooth  and  turbulated  sharp  180  degree  turns.  For 
gradual  180  degree  turns,  the  pressure  loss  coefficients 
were  based  on  the  equations  from  reference  [25]  and 
the  correlations  for  the  heat  transfer  coefficient  vari¬ 
ation  was  based  on  the  experimental  data  of  Baughn, 
et  al.  [20]. 

c)  Finned  Passages:  For  cooling  passage  channels 
with  fins  parallel  to  the  flow  direction  (usually  near 
the  trailing  edges  of  blades,  figure  1(b)), an  effective 
heat  transfer  coefficient  is  obtained  [3]  which  can  be 
applied  to  the  corresponding  plain  passage  surface 
area  without  fins. 

f)  Flow  Branchim;:  Pressure  loss  due  to  flow  branch¬ 
ing  has  been  modeled  based  on  reference  [25]. 

r)  Leading  Edge  Tin  Can  Impingement  Cooling: 
Modeling  of  the  tip  cap  impingement  cooling  is  based 
on  the  development  by  Meitner  [3]  who  utilized  the 
data  of  Chupp  et  al.  |27].  It  includes  the  effects  of 
rotation  on  pressure  clianges  across  the  tip  cap  im¬ 
pingement  chamber  and  the  temperature  increase  due 
to  impingement  in  the  chamber.  The  model  also  com¬ 
putes  the  mass  flow  through  the  impingement  holes. 
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Corrections  have  been  applied  to  the  Meitner  model 
based  on  the  work  of  Epstein  et  al.  [28]  which  include 
the  effects  of  rotation  on  heat  transfer  in  a  impinge¬ 
ment  cooled  turbine  blade.  Modeling  of  by-pass  flows 
(see  figure  1  (b))  is  also  based  on  Meitner’s  develop¬ 
ment  [3|. 

h)  Rotational  Effects  of  Plain  Passages:  Effects  on 
both  pressure  drop  and  heat  transfer  coefficient  due 
to  Coriolis  forces,  centrifugal  forces,  and  centripetal 
buoyancy  have  been  included  for  ducts  (circular, 
square  or  rectangular)  rotating  in  orthogonal  and/or 
parallel  modes  (with  respect  to  the  axis  of  rotation) 
and  for  both  laminar  and  turbulent  flows. 

For  laminar  flow  in  plain  passages  with  smooth 
or  rough  walls,  and  with  their  axes  normal  to  the  axis 
of  rotation,  the  pressure  loss  correlations  are  based  on 
the  work  of  Mori  and  Nakayama  [29]  and  heat  trans¬ 
fer  correlations  are  based  on  the  work  of  Morris  [30] 
and  for  moderate  temperature  differences  on  the  an¬ 
alytical  work  of  [29].  For  turbulent  flows,  the  heat 
transfer  correlations  cire  based  on  the  work  of  Clifford 
[31]  or  Morris  and  Ayhan  [32,33]  with  corrections  for 
the  leading  and  trailing  surfaces  based  on  Harsagama 
and  Morris  [34].  The  user  also  has  the  choice  of  bas¬ 
ing  it  on  the  theoretical  work  of  Mori,  Fukada  and 
Nakayama  [35],  The  pressure  loss  correlations  are  ob¬ 
tained  from  the  work  of  Ito  and  Nanbu  [36]  or  from 
reference  [35]. 

For  laminar  flow  in  plain  passages  rotating  about 
a  parallel  axis,  the  pressure  drop  correlations  are 
based  on  the  rotating  friction  factor  given  by 

frot  —  f»t  +  fez  (12) 

/„  =  0.503  Ro*  (13) 

where  is  “the  excess  friction  factor”  suggested  by 
Morris  [36],  The  heat  transfer  correlations  are  based 
on  the  work  of  Morris  and  Dias  [38]  and  for  moder¬ 
ate  temperature  differences  on  the  theoretical  work 
of  Mori  and  Nakayama  [39].  For  turbulent  flows,  the 
heat  transfer  correlations  are  based  on  the  work  of 
Morris  and  Dias  [40]  and  Woods  and  Morris  [41]  and 
for  moderate  temperature  differences  on  the  theoret¬ 
ical  work  of  Nakayama  [4^.  The  pressure  drop  corre¬ 
lations  are  based  on  [43].  For  the  entrance  region,  the 
heat  transfer  correlations  are  based  on  work  of  Morris 
and  Woods  ]43]. 

It  has  been  shown  from  the  experimental  data 
of  Hajek  et  al.  [44],  and  Guides  [45]  that  the  heat 
transfer  characteristics  due  to  rotation  are  radically 
different  for  the  leading  and  trailing  surfaces  and  also 
for  the  centrifugal  (radially  outward)  and  centripetal 
(radially  inward)  passages.  These  features  have  been 
captured  in  the  correlations  used  in  the  present  code. 

i)  Rotating  Channels  with  Heat  Transfer 

Enhancement  Devices: 


transfer  coefficient  and  friction  factor  (pressure 
drop)  distribution  due  to  enhancement  devices  and 
then  computing  the  effects  of  rotation  as  a  ra¬ 
tio  of  a)  Nu^otating)/Nu(non-rotating)  and  b) 
f(rotating)/f (non-rotating)  based  on  the  experimen¬ 
tal  correlations  mentioned  above.  The  validity  of  the 
correlations  used  is  shown  in  the  following  section. 

It  should  be  noted  that  the  user  always  has  the 
option  of  providing  the  pressure  loss  coefficient,  fric¬ 
tion  factor,  and/or  heat  transfer  coefficient  for  any 
given  interval  in  the  form  of  data  to  be  curve  fitted, 
equations,  or  a  set  of  given  values. 

COMPARISON  WITH 
EXPERIMENTAL  DATA 

In  the  following  paragraphs,  the  predictons  from 
the  code  are  compared  with  a  wide  variety  of  exper¬ 
imental  data  for  various  types  of  heat  transfer  en¬ 
hancement  devices.  As  mentioned  earlier,  the  predic¬ 
tions  from  the  code  are  compared  with  the  data  from 
the  sources  which  are  different  from  those  used  in  the 
correlation.  Both  non-rotating  and  rotating  passage 
data  with  and  without  various  heat  transfer  enhance¬ 
ment  devices  are  compared  to  give  a  clear  picture  of 
the  range  of  application  of  the  code. 

a)  Non-Rotating  Passage  Data: 

For  a  straight  duct  with  90  degree  turbulators, 
the  code  predictions  are  compared  with  the  experi¬ 
mental  data  of  Taslim  and  Spring  [46,47]  in  figures 
2  and  3.  Figure  2  compares  the  mean  Nusselt  num¬ 
ber  variation  with  Reynolds  number  (Re)  for  a  fixed 
rib  height  to  channel  hydraulic  diameter  ratio  (e/Dh, 
also  known  as  blockage  factor),  rib  pitch  to  height  ra¬ 
tio  (P/e)  and  channel  aspect  ratio  (W/H).  The  defini¬ 
tion  of  these  terms  are  clarified  in  the  inserts  in  figures 
2  and  3.  Figure  3  compares  the  Nusselt  number  en¬ 
hancement,  expressed  as  (Nu/7Vu,f),  with  a  variation 
in  the  blockage  factor  e/Dh  for  fixed  values  of  P/e, 
W/H,  and  Re.  Predictions  are  in  excellent  agreement 
with  the  experimental  data. 

In  figures  4  and  5,  the  code  predictions  for  a  two- 
pass  serpentine  passage  with  a  180  degree  turn,  for 
both  smooth  and  turbulated  walls,  are  compared  with 
experimental  measurements.  The  calculated  Nusselt 
number  variation  along  the  passage  is  compared  with 
the  data  of  Chandra  ]48]  in  figure  4  and  with  the 
data  of  Boyle  ]11]  in  figure  5.  Agreement  with  data 
of  ]48]  is  very  good  for  both  smooth  and  turbulated 
passages.  The  reason  for  the  reduction  in  the  Nusselt 
number  in  the  180  degree  turn  region  (X/D/,  =  12 
to  14)  for  the  turbulated  case  is  that  there  are  no 
turbulators  in  the  turn  region.  In  the  case  of  Boyle’s 
data,  the  180  degree  turn  region  (X/Dh  =  14  to  16) 
has  turbulators. 

Predictions  agree  very  well  with  Boyle’s  data  ]ll[ 
at  higher  values  of  Re  where  as  at  low  values  of  Re  (< 
=  20000),  the  agreement  is  good  except  in  the  region 


The  combined  effects  of  rotation  and  heat 
transfer  enhancement  devices  (like  ribs,  pin  fins 
etc.)  are  obtained  by  first  computing  the  heat 
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Figiirc  2.  Coiuj).iri?oii  of  Nussolt  number  variation 
with  Reynohls  number  witli  the  turbulated  passage 
data  of  Taslim  and  Spring  ('Ifi.'l?) 
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Figure  3.  Comparison  of  Nusselt  number  variation 
with  blorkage  fartor  [e/D),)  with  the  tubulated  pas¬ 
sage  (lata  of  'Faslim  and  Spring  146,47]) 


54  —  cc-o  ITjBBi  -  30000 

1  3  £*0  lTo«3J  *•'2  “  -  '^63 


0  5  10  15  20  25  30 


Non-dlmenslonol  distance  oleng  f!c"*  poth,  X/D 

h 


Figure  4.  Comparison  of  variation  of  Nusselt  number 
with  X/F);,  with  the  two-pass  rib  roughened  passage 

fat  90"  with  a  sharp  J80  deg.  turn)  data  of  Chandra 
48|. 
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Figure  5.  Comparison  of  variation  of  Nusselt  number 
witJi  X/ Uh  with  the  two-pass  turbulated  passage  data 
of  Boyle  [1  Ij. 

immediately  after  the  180  degree  turn  where  the  com¬ 
puted  Nu  drops  sharply  compared  to  the  exjierimen- 
tal  data.  'I'he  maximum  error  is  of  the  order  of  20  to 
TO%  in  this  region.  The  reason  for  this  behavior  is 
not  known  at  present. 

For  pa.ssages  with  pin  fins,  the  code  predictions 
are  compared  with  the  experimental  data  of  ban  et 
al.  [49]  and  Feng  |50j  in  figures  6  and  7.  Figure  0 
compares  the  predicted  Nu  variation  with  row  numher 
with  the  data  of  ban  [49]  and  figure  7  rompares  the 
predicted  variation  of  friction  factor  (f)  and  Colburn 
factor  (j)  with  neynolds  number  with  the  data  of  Fong 
|50|.  Both  the  figures  show  that  the  predictions  arc 
in  excellent  agreement  with  the  experimental  data. 
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Figure  8.  Sclieiiiatic  of  tlic  rotating  serpentine  pas¬ 
sage  of  Ilajek  (44]  showing  tlie  four  soncs  compared 
with  tlie  code  predictions, 


Figure  0.  Comparison  of  variation  of  Nuaselt  number 
witli  row  number  willi  the  pin  fln  data  of  Lau  (49|. 
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Figure  9(a)-(d)  compare  the  Nusselt  number 
variation  with  uossby  number  (Ro)  for  "the  above 
mentioned  four  zones  along  the  passage.  As  seen  in 
Ggures  9(a)  and  0(b),  the  predictions  clcarl)'  sliow  the 
experimentaliy  observed  trend  for  llie  centrifugal  sec¬ 
tions  tones  1  and  2,  that  as  the  Uossby  number  in¬ 
creases,  due  to  increasing  rotational  speed.  Nusselt 
number  increases  on  the  trailing  (presoure)  surface 
and  decreases  on  the  leading  (suction)  surface.  This 
trend  is  observed  for  both  smooth  and  turbulatcd  pas¬ 
sages,  This  is  due  to  the  increased  secondary  flow 
effects  at  higher  values  of  Ro.  Figures  9(c)  and  0(d) 
compare  Nusselt  number  variation  with  no  for  the 
centripetal  channel.  These  figures  also  show  that  the 
trend  of  Nu  variation  with  Ro  is  in  general  opposite  to 
that  of  the  centrifugal  section.  Comparison  of  pres¬ 
sure  variations  along  the  p.issage  with  the  data  of 
Hajek  |44j  is  made  in  Ggure  10  for  the  turbulatcd  pas¬ 
sages.  All  the  above  figures  clearly  show  that  the  code 
predictions  agree  well  with  the  data  of  Hajek  [44], 


Figme  "  Comparison  of  variation  of  friction  factor 
and  Coiowwi  lUciui  wall  Reynolds  number  with  the 
pin  fill  data  of  Peng  |50|, 

1>)  Rotating  P,%sfiage  Data:  Extensive  comparison  has 
been  made  with  the  rotating  passage  data  fbotli 
amootli  and  turbulated  serpentine  passages)  of  Hajek, 
et  al.  l44[.  A  schematic  of  the  serpentine  passage  of 
Hajek  is  shown  in  figure  8.  As  shown  in  Ggure  8,  four 
typical  zones  have  been  chosen  for  comparison  with 
code  predictions  to  illustrate  the  validity  of  the  corre¬ 
lations  developed  for  rotating  channels.  Zones  1  and 
2  arc  in  the  radially  outward  (centrifugal)  section  and 
zones  3  and  4  are  in  the  radially  inward  (centripetal) 
Bcctioii  of  the  passage.  As  explained  earlier,  separate 
correlations  are  used  to  predict  the  effects  of  rotation 
on  the  leading  surface  (suction  surface),  trailing  sur¬ 
face  (pressure  surface)  and  the  two  side  walls  of  the 
coolant  passage,  This  enables  the  code  to  compute 
the  rotational  effects  more  accurately. 
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(a):  Zone  1. 

Figure  9,  Comparison  of  Nusselt  number  variation 
with  Rossby  Number  for  the  smooth  and  45  deg. 
turbulated  serpentine  rotating  passage  data  of  Hajek 
(44|. 


7 


NJ/NL.55 
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Non-dlmen»lonol  <JI«topce  along  flow  path,  X/0^ 


Kigtiie  10.  CoinpariMO!)  of  preimtie  (1iKlril>iition  .ilong 
llie  p.iHs.nge  will)  tlie  tnrbiilalcd  rotaliiig  ecrpcnlinc 
p.-uteage  data  of  llajck  \M\. 

C/oin)>nriaon  of  Nneaelt  )iuiiil>er  variation  will)  the 
centrifugal  paesage  data  of  Gnidcz  |'iC]  ia  made  in  fig¬ 
ure  11  and  with  the  data  of  lakakov  mitl  'IViiahin  \f>l\ 
figure  12,  Thcae  fignrca  ahow  that  the  rode  prcdic- 
tiona  are  in  good  agreement  witli  the  cxporimcntal 
data  in  both  caaca  with  a  maximum  error  of  about 
20%,  The  Ircnda  ahown  in  theac  fignrea  are  the  aamc 
aa  the  centrifugal  paaaagc  data  of  llajek  \i4]. 
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Figure  11  ,CJnui|>nf Jeon  of  Niiw^ell  niiinbcr  varintioii 
wliJi  lloefll)y  Ntiiiiher  with  the  flnta  of  OiihlcR 
at  X//;,,  ^  7.<1. 


(d):  Zone  i. 
Figure  9,  Concluded 
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Figure  12.  Comnarison  of  Nueselt  number  witli 
tlic  parameter  Gr/Re  with  the  data  of  Iskakov  and 
IVushin  jSl]  at  X/D/,  =  7. 

COOLED  RADIAL  TURBINE  ROTOR  ANALYSIS 

The  coolant  flow  of  a  cooled  NASA  research  ra¬ 
dial  turbine  rotor,  which  will  be  tested  at  the  Lewis 
Research  Center,  has  been  analysed  using  the  code. 
The  cooling  passage  consists  of  a  main  passage  split¬ 
ting  into  three  branches  and  includes  pin  fins  and 
finned  passages.  The  schematic  of  the  cooling  pas¬ 
sage  is  shown  in  figure  1(c).  For  a  given  inlet  total 
temperature  and  pressure  and  exit  pressures  from  the 
branches,  the  juesent  code  was  used  to  predict  the  dis¬ 
tribution  of  pressure  along  the  entire  passage  and  also 
the  mass  flow  through  tlie  main  and  sub  branches. 
The  results  are  shown  in  figure  13  and  table  1  respec¬ 
tively.  These  predictions  will  be  compared  with  the 
experimental  results  as  they  become  available. 


Table  1.  Mass  flow  computatiojis  for  the  coolant  pas¬ 
sage  of  a  NASA  researcli  radial  turbine 


Branch 

Code 

Computation 

Main 

6,734 

B1 

2.399 

B2 

1.C32 

B3 

2.703 

riv,/  =  O.OOlkg/s 

WAIN  rJASSAOr. 


Flow  path  location 


Figure  13.  Pressure  distribution  along  the  coolant 
passage  of  a  NASA  research  radial  turbine. 

CONCLUDING  REMARKS 

A  generalized  one  dimensional  flow  code  appli¬ 
cable  to  compressible  flow  through  complex  rotat¬ 
ing  turbine  cooling  passages  with  various  heat  trans¬ 
fer  enhancement  devices  has  been  developed.  Even 
though  this  is  a  1-D  code,  the  effects  of  rotation  on 
Iieat  transfer  is  taken  into  account  in  a  quasi  2-D  fash¬ 
ion,  i.e.  a  separate  correlation  is  available  for  eacli 
side  of  a  rectangular  passage.  The  code  has  been  val¬ 
idated  using  a  wide  variety  of  both  non  rotating  an<l 
rotating  passage  data.  The  authors  feel  that  the  code 
will  be  useful  to  the  turbine  designer  in  evaluating 
the  effects  of  various  cooling  j)assage  designs  on  the 
overall  performance  of  the  turbine.  For  more  accurate 
compulations,  it  may  be  necessary  to  include  a  local¬ 
ized  3-D  analysis  where  flow  reversals  and  flow  stag¬ 
nation  may  cause  the  1-1)  code  predictions  at  these 
locations  to  be  in  error  by  a  wide  margin. 

NOMENCLATURE 
A  area 

Cp  specific  heat  at  constant  pressure 

D  diameter 

Df,  hydia\ilic  diameter  defined  by 

L>u  -  4A/P 
d  pin  fin  diameter 

e  surface  roughness;  or  rib  height 

f  fanning  friction  factor  defined  by 

ip Ih./{21V^P) 

Gr  rotational  Grashof  number  defined  by 
Gr  =  - 

H  flow  cliannel  height 

h  enthalpy 

h„  coolant  heat  transfer  coefficient 

j  Colburn  factor  defined  by 

j  -  h.rr^/y(rVCp) 
k  thermal  conductivity 

f.  inteival  length 

til  coolant  tiia.ss  flow  rate  at  inlet  to 

the  channel 
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N  rotational  speed 

Nu  Nusscit  number;  Nu  =  hel^D^k) 

NUft  Nusselt  number  for  a  stationary  plain 
passage 

P  perimeter;  or  pitch  of  turbulators 

Pr  Prandtl  number;  Pr  =  Cp^i/k 

p  pressure 

Prej  reference  pressure 

R  gas  constant 

Ra  Rotational  Rayleigh  number;  defined  by 
Ra  ^  GrPr 

-  -  TaJ)DlPTMp)'^ 

Rc  Reynolds  number;  defined  by 

Re  =  pVDjp 

Ro  Rossby  number;  defined  by  Ro  =  ND^/V 
T  radius;  co-ordinate  direction  used  in 

the  code 

S  pin  fin  spacing  in  spanwise  (transverse) 

direction 

T  temperature 

Tf,  bulk  temperature  of  the  fluid 

Tu,  wall  temperature 

Tnw  adiabatic  wall  temperature 

V  fluid  velocity 

W  flow  channel  width 

W/II  channel  aspect  ratio  (ribs  (if  any)  on 

side.  W) 

X  distance  along  the  flow  path 

X  pin  fin  spacing  in  streamwise  direction  or 
streamwisc  distance  from  entrance  to 
channel 

y  co-ordinate  direction  used  in  the  code 

z  co-ordinate  direction  used  in  the  code 

ft  volume  expansion  coefiTicient 

7  ratio  of  specific  heats 

5  a  change  in  a  quantitiy 

/X  absolute  viscosity 

•  p  density 

SUBSCRIPTS 

aw  adiabatic  wall 

cx  excess  factor  due  to  rotation 

inj  injection 

n  interval  number 

et  stationary  case 

ref  reference  value 

rot  rotating  case 

X  component  in  the  direction  of  main 

stream  flow 

SUPERSCRIPTS 

(  )’  total  conditions 
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